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Abstract 

We consider a discrete set of individual agents interacting with a continuum. Examples 
might be a predator facing a huge group of preys, or a few shepherd dogs driving a herd 
of shceps. Analytically, these situations can be described through a system of ordinary 
differential equations coupled with a scalar conservation law in several space dimensions. 
This paper provides a complete well posedness theory for the resulting Cauchy problem. 
A few applications are considered in detail and numerical integrations are provided. 
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1 Introduction 

In various situations a small set of individuals has to interact with a continuum. A first famous 
examples comes from the fairy tale of the pied piper [7|, where a musician frees a city from rats 
using his magic flute. An entirely different case is that of shepherd dogs confining sheeps while 
pasturing, or that of a wild predator seeking to split a flock of preys. From a deterministic 
point of view, studying these phenomena leads to a dynamical system consisting of ordinary 
differential equations for the evolution of the agents and partial differential equations for 
that of the continuum. Here, motivated by the present toy applications, we choose scalar 
conservation laws for the description of the continuum's evolution. In particular, no diffusion 
is here considered. On one side, this choice makes the analytical treatment technically more 
difficult, due to the possible singularities arising in the density that describes the continuum. 
On the other hand, we obtain a framework where all propagation speeds are finite. As a 
consequence, for instance, a continuum initially confined in a bounded region will remain in 
a (larger but) bounded region at any positive time. This allows to state problems concerning 
the support of the continuum, such as confinement problems (the rats should leave the city, 
or the shepherd dogs should keep sheeps inside a given area) or far more complex ones (how 
can a predator split the support of the density of its preys?). 

In the current literature, similar problems have been considered with a great variety of 
analytical tools, see for instance [2] for a fire confinement problem modeled through differential 
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inclusions, or |3j for a tumor-induced angiogenesis described through a stochastic geometric 
model. Other examples are provided by the interaction of a fluid (liquid or gas) with a rigid 
body or with an elastic structure, like a membrane, see [X2^ I13J : the evolution of the rigid 
body is described by a system of ordinary differential equations, while the evolution of the 
fluid is subject to partial differential equations like Navier-Stokes or Euler equations. Further 
results are currently available in the ID case. For instance, a problem motivated by traffic 
flow is considered in [9]; the piston problem, a blood circulation model and a supply chain 
model are considered in pQ. 

Formally, we are thus lead to the dynamical system 



where the unknowns are p and p. The former one, p = p(t, x) is the density describing 
the macroscopic state of the continuum while the latter, p = p(t), characterizes the state 
of the individuals. It can be for instance the vector of the individuals' positions or of the 
individuals' positions and speeds. The dynamics of the continuum is described by the flow /, 
which in general can be thought as the product / = pv of the density p and a suitable speed 
v = v(t,x, p,p). The vector field <p defines the dynamics of the individuals at time t and it 
depends from the continuum density p{t) through a suitable average A (p(t)) . Our driving 
example below is the convolution in the space variable A = p(t) * r], with a smooth 

compactly supported kernel r\. 

Below we address and solve the first mathematical questions that arise about (|f i.e. the 
existence and uniqueness of entropy solutions, their stability with respect the data and the 
equation, and the existence of optimal controls. A first well posedness result, that applies to 
general initial data, is provided in Theorem 12 .21 As usual in this context, see also [H l5| 18} [TT] . 
the hypotheses on / are rather intricate. However, the present framework naturally applies 
to situations in which the continuum can be supposed initially confined in a bounded region, 
i.e. p vanishes outside a compact subset of M. Nx . In this case, Corollary 12.51 applies and the 
hypotheses on / are greatly simplified. 

The next section presents the analytical well-posedness results. Section [3] is devoted to 
various applications, while all proofs are deferred to the last section. 

2 Notation and Analytical Results 

We now collect the various assumptions on that allow us to prove well posedness, i.e. the 
existence of solutions, their uniqueness and their stability with respect to data and equations. 
The hypotheses collected below are essentially those that ensure the well posedness of the 
conservation law and, separately, of the ordinary differential equation. 

Throughout, we denote R + = [0, +oo[ and B R N p (x,r) denotes the closed ball in ~R Np 
centered at x with radius r. Let T max G [0, +oo] and call / = [0, T max ] if T max < +oo, while / = 
R + otherwise. The real parameter R, i.e. the maximal possible density is fixed and positive. 
For a given compact set K in R v and a T > 0, we denote Qt = [0,7"] x R Nx x [0, R] x K. 




;) G R + x R Nx 
p G R + 
p G M N " 
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Flow of the continuum: at point x and time t, the continuum flows with a flux / = 
/ (t, x, p(t, x),p(t)} that depends on time t, on the space variable x, on the continuum density 
p evaluated at (t, x) and on the state p of the individuals at time t. We require the following 
regularity: 

(f) The flow f: I x R Nx x [0, R] x R n p -»• R N » is such that 
(f.l) / G C 2 (/ x R Nx x [0,R] x R^R^*). 

(f.2) For all (t,x,p) Elxl^'x R^p, f(t,x,0,p) = f(t,x,R,p) = 0. 

(f.3) For all T G / and for all compact subsets K C R p , there exists a constant 
such that for t G [0, T], s G R^, p G [0, i?] and p £ K, 

\\dpf(t,x,p,p)\\ < Cf, \div x f(t,x,p,p)\ <C f . 

(f.4) For all T G / and for all compact subsets -fT C R p , there exists a constant (7/ 
such that for t G [0, T], x G R^, p G [0, i2] and p G if, 

||V x .ap/(t,3;,/),p)|| < C/ . 
(f.5) For all compact subsets K C R Np , there exists a constant Cf such that 



sup \\V x div x f(t,x,p,p)\\dxdt < Cf , 

I JR N x p£K,p£[0,R] 

(f.6) For all compact subsets K C R p , there exists a constant such that 



sup \\div x f(t,x,p,p)\\ dxdt < Cf . 

(f.7) For all T G / and for all compact subsets -RT C R^ p , there exists a constant Cf 
such that for t G [0, T], p G [0, i?] and p£ K, 



/ ||V p div a; /(t,x,p,p)|| dx < C/ , ||Vp3p/(t,x,yO,p)|| < C/ for all a; G 



DAT- 



Condition (fj2]) states that at the maximal density p = R, the continuum is at congestion 
and can not move. Assumption (f J2J) has a key importance. The first part ensures the finite 
propagation speed of the solution to the partial differential equation, see Proposition 14. 21 or [8 , 
Theorem 1]. The second part ensures that the solutions are bounded, similarly to the role of 
the sublinearity {<p-3) in the ordinary differential equation. 

All these assumptions are satisfied, for instance, by vector fields of the form u(p,x,p) = 
v(p)v(x,p) with v G C 2 ([0, J R];R) and v G C 2 (R^ x R V *;R JV '*). 

We note that if / does not depend explicitly on t and x, which is a usual situation when 
dealing with systems of conservation laws in one space dimension, then the above assumptions 
reduce to only (fUJ), (fJ2D, the first part of (fJ3J) and the first part of (f^TJ). 

Moreover, Corollary 12.51 shows that whenever the initial density distribution p has com- 
pact support, then the requirements on / are reduced, since only (f^TJ), (fJ2D and (f J3J) are 
necessary. 
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Speed of the individuals: at time t, the individuals' state changes with a speed ip = 
(p (^t,p(t),A (p(i)) (p(t))^j that depends on time t, on the individuals' state p at time t and 

on an average A (p(i)) of the continuum density p evaluated at time t and computed at p(t). 
On the averaging operator A we require the following conditions. 

(A) A: L^R^M) ->■ W 1 - 00 ^^; R Nr ) is linear and continuous, i.e. there exists a constant 
C A such that for all p G L^R^jK) 

PpI| w i.°° ^ c a ||p|| L i • 

Below, the operator norm of A is denoted [|-A||^rLi w 1 * 00 )- For instance, in the case N p = N x , 
a typical example of such an operator A is (A(p)) (p) = (p*r])(p) for a kernel r\ G C^R 7 ^; R) 
with f RNp rjdx = 1. 

The speed law tp satisfies the assumptions: 

(<p) The vector field tp: R + x R^ x R^ R n p is such that 

(<£>.l) 1 1—?- ip(t,p,r) is measurable for all p G R^ and all r G R^; 

(tp.2) there exists a function C v G L 1 (/;M + ) such that for a.e. i G /, Pi,P2 £ K p and 
n,r 2 G R Nr , 

||v?(t,pi,ri) - ^(t,P2,r 2 )|| < C ¥ {t) (||pi -P2II + ||n - r 2 ||) ; 

(<£>.3) there exists a function G L 1 (/;R + ) such that for a.e. t G [0,T], for all p G R"^ 
and for all r G 1* , 

\\v(t,P,r)\\ < C v (t) (1 + IHI) . 

These hypotheses are motivated by the standard theory of Caratheodory ordinary differential 
equations, see [3 § 1]. All the above assumptions (f), (A) and (cp) are satisfied in the 
applications considered in Section O 

As a first step in the analytical treatment of we rigorously state what we mean by 

solution to (jl.ip . 

Definition 2.1 Fix p G (L 1 n BV) (r n *; [0,R]^j and p G R^V A pair (p,p) with 

p G C° (^I; 'L 1 (R Nx ; [0, R])^ and p e W 1 ' 1 (I;R Np ) 

is a solution to with initial datum (p,p) if 

(i) the map p = p(t,x) is a Kruzkov solution to the scalar conservation law 

d tP + div x f(t,x,p,p(t)) =0 (2.1) 

(ii) the map p = p{t) is a Caratheodory solution to the ordinary differential equation 

p = <p(t,p,A(p(t))(p)y, (2.2) 

(Hi) p(0) = p and p(0) = p. 
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For the standard definition of Kruzkov solution we refer to [H Definition 1], for that of 
Caratheodory solution, see [SI § 1]. 



Theorem 2.2 Under conditions (f), ((p) and (A), for any initial datum p £ R p and p £ 
(L 1 nBV)(M Ar -; [0,i?]), problem ^JD admits a unique solution in the sense of Definition \2.1\ 
This solution can be extended to all I . 

Let now f\, f% satisfy (f); A\, A<i satisfy (A) and <pi, if2 satisfy (<f); in all cases for the 
same interval I and the same parameters or functions R,Cf,CA,C v . Then, for any initial 
data (pi,pi),(p2,ih) G (L 1 nBV)^ 1 ; [0, R]) x R n p , the solutions (pi,Pi) and(p 2 ,P2) to the 
problems 



d t pi + div x fx (t, x, pi , pi(t)) =0 
Pi = fi (t,pi, (Aipx(t)) (pi) 
pi(0,x) = pi(x) 
Pl(0) =Pi 

satisfy the inequalities 

\\(pl -p2)(t)|| L i 

< (i + JC(t)) \\pi-fhWv 

+m (||W-/a)||L- (0t ) 



and 



< 



and 



d t p2 + div x f 2 (t, x, p 2 ,P2(t)) = 

P2 = ¥2 (t,P2, (A 2 p2(t)) (>2) 

p 2 {0,x) = p 2 {x) 

P2(0) =P2 



(2.3) 



+ ||div(/i-/ 2 )| 



L 1 



x)x~L°°(lO,t]x[0,R]xKt) / 



+JC(t) [\\ipi - Hll/*>([0,t]x.fir t x[0,C' A ]) + Pi ~ ^H^L^wi.oo) + ||pi 



P2\ 



\\(Pl -P2)(t)\\ 

(l + £(t)) ||pi-p 2 || 

+m(\\o P (h-f2)\\^ {nt) 



+ ||div(/a-/ 2 



x)x~L°°([0,t]x[0,R]xK t ) 



+JC(t) [\\<pi - y2|| Lo0 ([o )t ]x^ t x[o,CAD + W Al ~ A2 Wc(L\wi>°°) + Wpi - y°2 || L i 
where K, £ C°(/;R + ) vanishes att = 0. 

More detailed expressions of the various coefficients are presented in Section [H The proof, 
which is deferred to Section 21 is obtained through Banach Contraction Theorem. The nec- 
essary estimates for the convergence are a consequence of [81 Theorem 5], [5l Theorem 2.5] 
and of an adaptation of the standard theory of Caratheodory ordinary differential equations, 
collected in the following two lemmas. 



Lemma 2.3 Let (f) hold. Choose any p £ (L 1 n L°° n BV)(ffi w - [0,R]). Fix a function 
7r 6 C°(I;R iVp ). Then, the conservation law 



d t p + div x f (t, x, p, 7r(*)) = 
p(0,x) = p(x) 



(2.4) 



admits a unique solution p G C° yL] L^R^, [0, R])j ■ For all t £ L, introduce the compact 
set K t = B k n p (0, ll 7r llc°([o t]))' denote Q t = [0, t] x R x x [0,R] x K t and define 

Kt = (2N x + l)\\V x d p f\\ Lao{nt) . (2.5) 
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Then, the following BV estimate holds: for all t G I 



TV (/>(*)) < (TV(p) + N a> W Nm tJ ll V - div -/(-^'-)|| LO c a0i , ]x[0ii?Jx ^ ) dx) e^. (2.6) 

£e£ now, for i = 1,2, pi be the solution to U.l\) corresponding to the initial datum pi and to 
the equation defined by 7Tj G C°(I;W Np ) and by fi, satisfying (f). Then, 



(pi -P2)(*)|| L i < IIpi-P2|| l i 



+tC(t) 



Kl - 7T2|| L o O([0it]) + \\d p (fi ~ /2)|| L oc (f7t) 
+ ||div (/l - /2)|| L i( R JV a; ) xIjOO ([Q )t ] x [0 ) JJ| x J!f t ) 



where C(t) depends on TV(pi), || V x <9 p /i || Loo(nt) , || V* div^ /i|| L i(K^)xL°°([o,t]x[o,fl]x^ t ) anrf 
||Vp3 p /2|| Loo(nt) , ||div x Vp/ 2 || L i( R ]v ;c ) xLO o([ 0)t ] x [ 0i/? ] xXt ), 

An explicit expression of C(t) is provided in (|4,ip . 

The estimates related to the ordinary differential equation are provided by the following 
lemma. 

Lemma 2.4 Let (cp) and (A) hold. Choose an initial datum p G W Np and fix a function 
r G C° (V; L 1 (R^ ; [0,^])") . Then, the ordinary differential equation 

(2.8) 



p = (p[t,p,A(r(t)) (p) 
p(0) = p. 



admits a unique solution p G W^^° (/; ) . T/ie following bound holds: 

\\p(t)\\ < (\\p\\ + l)eti C ^ dT -1. (2.9) 

Given two initial conditions pi,p 2 G R"^, too functions r\,r 2 G C° fl; L 1 (R-^ ; [0,i?])^, too 
speed /aws (pi,<p% satisfying (<p) and two averaging operators A\, A2 satisfying (A), define 

F(t) = (l + CA[|ri|| L co ([ o, t ] ;L i)) [ C ^ dT ■ ( 2 - 10 ) 

Then, 

\\(Vl-P2)(t)\\ 

< e F ^\\p 1 -p 2 \\ + j\ F ^- F ^\\Mr)-Mr)\\ Loa dr (2 _ n) 
+ jy®- F VC v (T) (c A \\( n -r 2 )(r)|| L1 + \\A 1 -A 2 \\ c{L1 ^ 1 , ao) \\r 2 (r)\\ hl ) dr. 

In the applications below, the support of the initial data is compact. Thanks to the 
finite propagation speed typical of conservation laws, this allows a major simplification in the 
assumptions of Theorem 12.21 



Corollary 2.5 Consider problem (jl.ip with f satisfying f/OP, (f$2$ and (f$3\). Let A sat- 
isfy (A) and ip satisfy (<p). If p vanishes outside a compact set, then problem (jl.ip admits 
a unique solution in the sense of Definition \2.l\ This solution can be extended to all of I. 
Moreover, the stability estimates of Theorem \2.2\ apply, provided both p\ and p 2 vanish outside 
a compact set. 
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3 Applications 

This section is devoted to a few sample applications of Theorem 12.21 While the unknown p 
keeps throughout the meaning of a scalar density, the state p of the individuals is the position 
of a single agent in § 13.11 it is a vector of several positions in § 13.21 and it becomes a 4-vector 
position-speed in § 13.31 

Numerical integrations are also provided in order to show the qualitative behavior of the 
solutions. In all cases, the Lax-Priedrichs method, see \10\ § 12.5], with dimensional splitting 
was used for the conservation law and Euler polygonals to integrate the ordinary differential 
equation. 

3.1 The Pied Piper 

As a first toy application we consider the situation described in [7J n. 246]. To lure rats away, 
the city of Hamelin (now Hamel) hires a rat-catcher who, playing his magic pipe, attracts all 
mice out of the city. In this case, p = p(t, x) is the mice density and p = p(t) is the position 
of the piper. Rats move with a speed v(p) v(p — x), with the scalar v and the vector v having 
the qualitative behavior in Figure [TJ More precisely, at density mice have the fastest speed 



V i 


P -> v(p) 


\\v\\ 


^^^^^ 


p q(p) 





R p 





\\x\\ 


R p 



Figure 1: Left, v is assumed C 2 and decreasing. Center, v describes the attraction felt by the mice 
towards the piper. Right, q accounts for the acceleration of the piper when surrounded by a high mice 
density. 

while at density R their speed vanishes. The term v accounts for the attraction of the mice 
towards the piper. The magic musician has a speed q(p * rj)ip(t), i.e. he moves faster when 
the average density of mice around him is higher. On the contrary, when only few rats are 
near to him, he slows down. 

Lemma 3.1 Let N x = 2, N p = 2, N r = 1 and fix a positive R. Assume v E C 2 ([0, R]; R), 
v G C 2 (E 2 ;M 2 ) ; q G W 1 ' oo ([0, R}; M), tp G W 1 '°°(E+; R 2 ), rj G C2(M 2 ,]R) with J^rjdx = 1. 
Assume that v(R) = 0. Define 

f(t,x,p,p)=pv(p)v(p-x) <p(t,p,r) = g(r)^(t) Ap = p* x n. (3.1) 

Then, this setting fits in the framework of Corollary \2. 5\ as soon as p vanishes outside a 
compact set. 

The proof is immediate and, hence, omitted. 
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Numerical example: To fix a specific situation, we choose the following functions 
in (HU): 

v(p) = V max (1 - F max = 9 i? = 1 



<?(r) 
#) 

rj{x) 



x e 



COS Lot 

- sin Lot 



_3 

7rr. 



^ ^max |o, r p 2 — ||x|| 2 |^ 



w = 1 



r p = 0.15 



v p = 1 



(3.2) 



At time £ = 0, we assume that rats are uniformly distributed with density R = 1 in the 
rectangle [—0.5,0] x [0.35,0.85]. The piper starts moving at the point (—1,0.5). 




Figure 2: The pied piper and the rats, at times 0, when p = (0,0.5); 0.171, 0.543, 0.945, 1.447 and 
1.930, when the rats almost completely left the rectangle and p = (0.366, —0.983). 

Several optimization problems can now be stated with reference to fjl. lj) — f|3. lj) — (|3.2j) . Re- 
ferring to the situation [7, n. 246], a first natural question is the following. Let the compact 
set K be the area of the city and fix a finite positive time T max . Then, find the initial position 
p and the trajectory x]j of the piper so that the amount of mice left in the city at time T max 
is minimal. In other words, we want to minimize the functional 

(P, i>) J (yp{P, $)) (^max, X) dx . 

Here, p(p, i/>) is the p-component of the solution to (|l.ip - (|3.ip - (|3.2|) . The existence of such 
an optimal strategy for the piper follows from Theorem 12.21 via a standard application of 
Weierstrafi Theorem. 
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Proposition 3.2 Let T max be finite. Denote by K C R 2 the compact Hamelin urban area. 
Define the set of the possible piper's route choices 



K. = I (p, V) G K x W 1,00 (J; R ) : 



< 1 

■\Yl,oo 



and caZ/ : /C i— >■ R i/ie functional giving the total amount of mice in Hamelin at time 2~max; 
i.e. 

3{vA) = J (p(p~,$)) (T max ,x)dx , 

where p(p, rj)) is the solution to $l.l\) - $3. l\) - $3.2\) . Then, there exists an optimal trajectory 
(p*,V>*) G /C suc/z f/iai J"(p*,V>*) = minx; J" (p,^). 

Thanks to the stability estimates in Theorem l2.2| the proof of this proposition directly follows 
from Ascoli-Arzela Theorem that allows to prove the compactness of JC. 



3.2 Shepherd Dogs 

On the plane, consider a herd of, say, sheeps controlled by n shepherd dogs. Then, p is the 
density of sheeps and p = (pi, . . . ,p n ) is the vector of the positions of the dogs, so that each pi 
is in R 2 . We assume that initially the sheeps are distributed around, say, the origin and tend 
to disperse moving radially with a speed directed by v* r (x). The duty of the dogs is to prevent 
this dispersion and they pursue this goal moving around sheeps or, more precisely, with a speed 
(p orthogonal to the gradient of the sheeps' density. The sheeps modify their speed escaping 
from the dogs with a repulsive speed Vd(x,p) = v(x — pi), where v behaves qualitatively 
as in Figure [TJ Finally, the speed of the sheeps is then given by v(p) {v r {x) + ^22=1 v(x — Pi)) 
where v is maximal at the density zero and vanishes at the maximal density R. This last fact 
means that the sheeps can not move when their density is maximal. 



Lemma 3.3 Let n £ N ; N x = 2, N p = 2n, N r = 2n and fix a positive R. Assume v £ 
C 2 ([0,R};R), v r G C 2 (R 2 ;R 2 ) v G C 2 (R 2 ;R 2 ), n G C 2 (M 2 ,R) with J R2 r/dx = 1. Assume 
that v(R) = 0. Define 

f(t,x,p,p) = pv{p) (v r {x) + Yn=iv{x - Pi)) , 

^' r ) = yd 7iW' (3 - 3) 

Ap = p* x Vn. 

Then, this setting fits in the framework of Corollary \2.5\ as soon as p vanishes outside a 
compact set. 



Here, r = (r\, . . . , r n ) is a vector in (R ) n and we set r = 



l > 



with 



a 




b 


b 




—a 
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Numerical Integration: To fix a specific situation, we choose n = 2 and the following 
functions in (|1.1|) : 




v d = 100 

At time zero, sheeps are uniformly distributed at the maximal density R = 1 in the circumfer- 
ence centered at (0,0) with radius 0.2. Dogs start moving from (0.7, 0) and (—0.7, 0) Graphs 



• 

























Figure 3: Solution to CE~J]) "(El at times t = 0, t = 0.044, t = 0.067, t = 0.111, t = 0.156, 
t = 0.200. Sheeps are initially uniformly distributed at the maximal density R = 1 in the circumference 
centered at (0,0) with radius 0.2. Dogs start moving from (0.7, 0) and (—0.7, 0), they succeed in 
confining the dispersion of the sheeps, at least for the tie interval considered. 

of the corresponding solution are in Figure 

Merely technical modifications may allow to pass to various other problems. For instance, 
dogs may be asked to constrain the movement of all sheeps towards a certain area. 

3.3 Predator and Preys 

We consider here a predator attacking a group of preys. We think for example at a hawk 
pursuing a flock of smaller birds or at a shark attacking a group of sardines. Here, p is the 
density of the preys with x S M 3 , p is now the pair (P, V) where P S M 3 is the position 

of the predator, V G M 3 is its speed and we postulate below an equation for the acceleration 
P = V of the predator. Indeed, the framework in Theorem 12.21 allows to consider also second, 
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or higher, order ordinary differential equations for the single agents. The initial density of 
the preys is assumed to have a compact, connected support. The aim of the predator is to 
divide this connected group into two smaller (disconnected) groups. Hence, its acceleration 
is directed along the gradient of the preys' density, say P = ap(t) * x V77 for a suitable a > 0. 
The preys have a speed Vm ax (l — p/R)Vq, for a fixed Vq S M 2 , and modify it trying to escape 
from the predator. The resulting speed of the preys is thus 

v(t,x,p,p) = y max (l - p/R) (y + Be- C W x ~ p VW {x-p(t))) (3.5) 

where B,C are positive constants. The former one is related to the speed at which preys 
escape the predator and the latter to the distance at which preys feel the presence of the 
predator. Once again, v is maximal at zero density and vanishes at the maximal density R, 
which means that the preys can not move when their density is maximal. 



Lemma 3.4 Let N x = 3, N p = 6, N r = 3 and fix a positive R. Assume v is as in (|3.5p . 
r/ G Cc(M 2 ,M) with J R2 ndx = I. Denote p = (P,V) and define 

Ap = p* x Vr]. (3.6) 

Then, this setting fits in the framework of Corollary \2.5\ as soon as p vanishes outside a 
compact set. 



f(t,x,p,p) = pv(t,x,p,p), if it, 



p 

V 




Figure 4: Solution obtained through the numerical integration of (|1.1I) - (|3.5I) ~ (|3.6[) - (|3.8[) computed at 
times 0, 0.091, 0.267, 0.358, 0.449 and 0.491. Note that the predator succeeds in splitting the support 
of the preys. 
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Numerical Integration: For graphical purposes, we limit the numerical integration to the 
2D case. With reference to (|l.ip - (|3.5|) ~ (|3.6|) . we choose the following parameters 



y max = 2 , C = 5.25 , V Q = [0 - 0.5] T , B = 40 , A = 400 

o / [„ o „2l \ n c 

r„ = 0.5 . 



and the initial datum 

Pn = 



'})" 



(3.7) 





-0.8 



, V 



. ^( X ^) = X [ _0.2,0.2] (x)X [-0.2,-0.1]^ 



(3.8) 



The result is in Figure 01 Note that the predator succeeds in splitting the support of the 
preys. 



4 Technical Details 

Throughout, we let Wn x = J c , (cos 9) Nx d9. We state below a Gronwall-type lemma for 
later use. 

Lemma 4.1 Let the functions a G C°(L;R), p G W 1 ' 1 ^;^), 7 G C°(L;R+), A G C°(/;K) 
be such that 

rt 



A(t)<a(t)h(t) + J 7(r)A(r) dr 



Then, for all t E I, 
A(t) < a(t) 



/3(0)exp[ / a(r) 7 (T)dr j + ^ /3'(T)exp f y a(s) 7(5) ds j dr 



Proof. Using the following straightforward computations, we have: 

7(*)A(t) < a{t)P{t)j{t)+a(t)j(t) [ 7 (r)A(r)dr J 

JO 

[ e -J?«(ThWdT /" 7 ( T )A(r)dr ] < a(t) /3(t) 7 (t) e" ^ . 



\ jo 
Then, by integration we obtain 



/ 7(r)A( 
J 

Consequently, we have 



r ) d r < / a(t) p(t) 7(t) «M7«d« dr . 




A(t) < a(t) 



£(*) + / e^ a(s)7(s)ds a(r)/3(r) 7 (r)dr 



Integrating by part the last integral, we have finally the desired estimate. 



□ 
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Proof of Lemma \2.3l This proof consists in applying to the scalar conservation law 
dtp + div x f*(t, x, p) = with flux f*(t,x,p) = f (t, x, p,p(t)} first the classical Kruzkov 
result [U Theorem 5] and then the stability estimates in [5]. 

To apply Kruzkov Theorem, it is sufficient to verify condition (HI) in [5, Theorem 2.5] 
or the slightly weakened form in [llj. Note that: /* is C 2 in x and p by (f JT]) , and is C° 
in t by the regularity of it. This regularity is sufficient in the proof of Theorem 2.5], see 
also [8j Remark 4 in § 5] and [TT]. Moreover, for any t £ I 



d p f* G L°°([0,t] x R Nx x [0,i?];R^) and div x /* G L°°([0,i] x M JV * x [0,12]; 



Kruzkov Theorem can then be applied on any interval [0,t]. 

Observe that by (fj2]), the constant functions p(t,x) = and p(t,x) = R solve the con- 
servation law (|2.4|) . independently from tt. Then, by the Maximum Principle [5] Theorem 3], 
we have that any solution p to fj2.4|) satisfies p(t, x) G [0, R] for a.e. (t, x) G 1 x M. x and for 
all 7r G C°(l; R^). 

To prove the L 1 continuity in time and the TV bound, we apply [5> Theorem 2.5] in the 
weaker form in [TT]. To this aim, we verify also (H2) therein on the time interval [0,t], for 
any tel. By (fSJ) and the continuity of vr, V x d p f* G L°°([0,t] x R N * x [0, 12]; R N * xN *). 
Note also that, by (fj5]), Jq J^n x \\ Vx div x /*(r, x, p) || Loo dx dr < +oo, with an upper bound 
that depends on tt. 



We denote below = [0, t] x ] 



t>N x 



[0, 1?] x ETj where K t is as above. By [TT] Theorem 2.2] 



or [S] Theorem 2.5] we obtain the estimate 



TV (p(t)) < TV(p)e Ktt + N x W Nx f 

Jo 



/ V x div x f (t,x,-,7t(t)) 
. This implies ()2.6p . 



dx dr 



L°°([0,i?]) 

where k 4 = (2^ + 1)|| V a .5 p /|| LOO(Ct) . 

The L 1 -continuity in time of p follows from [5, Remark 2.4], thanks to (f]6]) and to the 
bound on the total variation, see also [U Proof of Lemma 5.3]. 

To estimate the dependence of the solution from the initial datum, we check the hypothe- 
ses (H3) in |TT] or [5] and apply [UJ Theorem 2.3] or [5l Theorem 2.6]. 

Let fi and f 2 satisfy (fJTJ), . . ., (fjS]). Assume that m and vr 2 are in C°([0, t], R n p). Let 
/* and /I be the corresponding compositions. With obvious notation, define K = K\ U 
and compute 

dpfi (t, x, p, tti(t)) - dpf % (r, x, p, 7t 2 (t)) 



< 



sup 

r£[0,t],a;eR JV ' ! < ! ,p6[0,ii] 

|9p/i - d ph\\ L oo (nt) + ||^V p / 2 || LOO(nt) ||tti - 7r 2 || L oo([o, t ]) 



which is bounded by (f]3]) and (fJTJ) . 

To complete (H3), it remains only to estimate the quantity 



< 



/7 

Jo Jw 
Jo Jw 



+ 



JJSL N x 



diva; (^/i (r,x,-,7Ti(r)) - / 2 (r, x, •, tt 2 (t) 

divsC/i - /a) (r,x,-,7n(r)) Loo dx dr 
V p div x / 2 (x)|| LOO ||7ri(r) - 7t 2 (t) || dxdr 



dx dr 



L°°([0,R]; 
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which is bounded thanks to (f j6|) and (f JTJ) . Now, we compare p\ and p2, obtaining 



|| (Pi -P3)(*)|| L i 
<I|P1 - P2|| L i 

+ 



TV(p) + N X W N . 

«i ' ./0 



* e «t(t— t) _ ]_ 



Va; div x . fi (r,x,-,vri(r)) 



L°°([0,fl]) 



dx dr 



x [\\dpfi - 9p/a|| L oo (nt ) + Fp V p/2|| Lt3 c (Ct) IKi - T2|| L <»([o,t]) 



+ 



^ div(/i - / 2 ) (r,a;,-,7ri(r)) 



/o Jm N x 

+ || V p d^ f 2 (r,x, 

which gives the final estimate. 



lL°°([0,,R]xi<- t ) 



L°°([0,fl]) 

7Ti(r) — (x) 1 1 ) dxdT 



(4.1) 

□ 



Proof of Lemma \2.4[ By ((f), we may apply [U theorems 1 and 2, Chap. 1] to (|2.8p and 
get the local in time existence and uniqueness of the solution. The bound (|2.9p follows from 
a standard application of Gronwall Lemma and ensures that the solution can be extended to 
the whole interval /. Assume for simplicity that ip\ and cp% satisfy (if) with the same function 
Cu>. Using the representation formula 



Pi=Pi+ I tpi(r,pi(T),A(ri(T)) (pi(r)fj dr 



we get 



\\(Pl ~ P2)(t)\\ 

< \\Pl -P2|| + / 
J 



<P1 [r,Pi(r),Ax (ri(r)) (pi(r))) - 922 (t,p 2 (t),^2 (?" 2 (t)) fe(r)) 



dr 



< WP1-P2W+ \\(<pi ~ y2)(r,pi(T),Ai(ri(r))(pi(T))) dr 

JO 

rt 



+ / ^ (r) _;p2)(t) ii + Ik 1 (ri(r) ^ (pi(r)) ■ ^ (r2(r) ) (p2(r)) 

< WPl-ihW+f C^(^)(l + ||V^i(ri)|| Loo ) ||(R-p 2 )(r)||dr 

£ ( L i jW i,oo)||(ri - ^2)(t)|| l i + 11^1 - yl 2||£( L i i w 1 .°°)|| r '2( T )||Li J ^ T 



\(<Pl ~ ^2)(t,-,-)\\raodT . 



An application of Lemma 14.11 with 



A(t) =||Pi-Pa||, 
a(t) =1 , 
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P(t) =\\p! - p 2 \\ + / \\(tpi - <P2)(t, •,•) Loo dr 

JO 

7(0 =C v (t) (l + Pi||£(Li,wi.-)ll r i|| L i 

+ / ll Al ll £ C L1 > W1 ' ao )IK ri ~ ^2)(t-) || L i + ||Ai -^2|| £ (Li,Wi>-)|| r 2(T)|| Ll dr. 

completes the proof of (|2.1ip . □ 
Proof of Theorem \2.2l The proof is divided in several steps. 

1. Local Existence. Here we rely on an application of Banach Fixed Point Theorem. Fix 
first the initial data p G (L 1 n BV)^; [0, R]) and p G R Np . Choose a positive f G I and, 
motivated by (j2.9|) . call 

J= (||p|| + l) e /o T ^W d ^-i. 

For any positive TZ, with J KJVa . pdx < TZ, and for any T G 0, T , define the complete metric 
spaces and the distance 



X p = \p£ C° ([CT]^ 1 ^;^^)) : sup I p(t, x) dx < TZ \ , 
{ te[o,T] Jr 

X = X p xC°([0,T};B RNp (0,5)) , 

d((pl,Pl);{p2,P2)) = SUp \\pi(t) - P2(*)|| L l + SUp \\pi(t) -p 2 (t)\\ 

te[o,T] te[o,T] 



Define the map T: X — > X by T(r, ir) = (p,p) if and only if p and p solve the problems 

(4.2) 



d t p + div x f {t,x,p,ir(t)) =0 ^ | p = ip(t,p,{Ar{t))(jp)) 
p(0,x) = p(x) <<U<1 J p ( ) =p. 



Note that both problems admit a unique solution, by lemmas 12.31 and 12.41 Moreover, by the 
conservative form of the former problem in (|4.2p . J RNx p(t, x) dx = f RNx p(x) dx < TZ, so that 
T is well defined. Moreover, Lemma [2.41 shows that the solution to the latter problem in (|4.2p 
is in W 1 ' 00 ([0,11; B RNp (0,5)) c C° ([0, T\; B RNp (0, S)) . 

To prove that T is a contraction, fix (fij7Ti) and (j2^2) and call (pi,Pi) = T(rj,7Tj). 
Then, define Kf = B r n p (0,5) and apply Lemma [2.31 with t = T. Note that Kt Q Kf. 

The former problem in (g^J) is then solvable in C° f [0, T]; L 1 (R iVa! ; [0, and the stability 

estimate (j2.7|) yields 

SUp ||pi(t) - P2(*)|| L i < TC(T) SUp 1 1 7T1 (* ) - 7T 2 (t)|| • 

te[o,T] J te[o,T] 

Apply now (|2.1ip 

sup ||pi(t)-j»2(t)|| < CU [ T C v (r)e F ^- F ^dT sup ||n(i) - r 2 (t)|| Ll , 

te[o,T] Jo *e[o,T] 

where i 7 is defined as in ()2.10p and can here be bounded as 

F(t) < (1 + C A ^) T C^(r) dr . (4.3) 
JO 
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Hence, 

d(T(pi,pi),T(p2,P2)) <max{rC(f),C A ( e F ( T )-l)} d {( Pl , Pl ), (p 2 ,p 2 )) . 

Choose now a sufficiently small T so that T is a contraction. Then, its unique fixed point is 
the unique solution to (jl.ip defined on the time interval [0, T]. 

2. Global Uniqueness: Let now (pi,pi) and (p 2 ,p 2 ) be two solutions to the same prob- 
lem (jl.ip and defined at least on a common time interval [0, T] C /. Define 

r* = su P {re [o,f] : (pi,pi)(i) = (p 2 ,p 2 )(t) for aine [o,r]} . 

By the uniqueness of the fixed point, (pi,Pi)(t) = (p 2 ,p 2 )(t) for all t E [0, T], so that the 
set in the right hand side above is not empty. Repeat Step 1 with initial datum (p*,p*) = 
(pi,pi)(T*) = (p2,P2)(T*), which is possible since p is bounded on [0, T*] and TV(p*) is 
bounded, by (|2.6p . Thus, we obtain that (pi,pi)(t) = (p 2 ,p 2 )(t) also on a right neighborhood 
of T* . This contradicts the maximality of T*, unless T* = T. 

3. Global Existence: Define now 



T* = sup {T E /: 3 a solution to (jl.ip defined on [0,T]} 
and assume that T* < +oo. By (|2.9|) . p is bounded on [0, T*[ and since 



||K*a) — K*x)|| < 



[ 2 C^{t) (l + \\p(t)\\) dr < [1+ sup [ 2 C^(T)dT 
Jti \ ■ / y te[o,T»] / Jti 



p is also uniformly continuous. Hence the limit p* = lirm T - p(t) exists and is finite. 

Apply now Lemma 12.31 on the interval [0,^], obtaining that the solution p to fj2.4j) is 
defined on all [0, T*] and, together with p, also solves (jl.ip . Now, we repeat Step 1 with 
initial datum (p*,p*) = (p, p) (T* ) , which is possible thanks to (|2.6p . In turn, this allows to 
extend (p,p) to a right neighborhood of T*. This contradicts the maximality of T*, unless 
T — T 

± * — ± max • 

4. Stability Estimates: Fix t > and let r E [0, t]. Let 1Z > max | J R]V;1 . p\ dx , J RjV:! , p 2 dx|. 
Then, by ([577|> and (f2TTj) . the solutions to (|iOj) satisfy 



||(pi -p2)(t)|| Ll 

< lIPl - P2|| L i + 



+ ||0o(/l-/ 2 ) 



+ ||div (/i - /2)|| L i( R jv ir ) xL ooQ 0)t ] x [o ) _R] x _f ! : t ) > 

||(Pl ~P2)(t)|| 

< e F ^\\p 1 -p2\\+ re F W- F M||(^-(p 2 )(r,-,-)|| T00 dr 



+ ^'e F(t) - F(r) C v (r) (c A ||(pi -p 2 )(r)|| Ll + ftpi - A 2 || £(L1)Wll00) ) dr. 
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with C as in LemmaESJ F as in (pL"gj) . if t = 5(0, 5 t ) and <J t = + l) e^o c ^ (r)dT - 1. Insert 
now the former estimate in the latter one and apply Lemma 14.11 with 

A = \\(Pi ~ P2)(t)\\ , 
a{t) = e F ®, 



+ / IK^-^^^OH^e-^dr + C^ f e- F ^C^T)\\- pi -p 2 \\^dT 
Jo Jo 

+C A f : rC(r)C,(r)e- f W 
Jo 

X - /2)|| L oo( Cr ) + ||dlV (/i - /2)|| L i( K JV a; ) ><L oo([ 0iT ] x [ 0>i? j x ^ T )) , 

7 (t) = C A t^(t)C(t)e F W, 
obtaining, with H(t, t) = exp J** C<p{s) (l + Ca 7£ + Ca sC(s)) ds, 

lbi-K|| < fexp ^F(i) + C a J*tC v <j)C(t) drjj ||p!-p 2 1| 

+ \J q n(T,t)drj hi ~ ^2\\ h ^ m]xKtX[0tCA ] 
+ \T^J o C v {T)n(T,t)dr \ \\Ax-A 2 \\ c(LljWl>00) 

+ [c A J*C lp (T)H(T,t)dTj ||pi-p 2 || Ll 

+ [c A J\c^r)c{T)n{T,t)d^ 

X JPp(/l - /2)|| L oo([ 0ji? ] xR JV a;xii - t ) + ||diVa, (/l - /2)|| L i( K iV a; ) ><L oo([ 0ii j] >< ^) 

Then, we immediately get the other bound 
llPi - P2|| L i 

< \\pi ~ fh\\ ^1 + tC{t) exp (^F(t) + C A £ tC v (t)C(t) dr 

+ (ntC(t) C v (T)«(T,*)drJ Pi -^2|| £(L i iW i^) 
+CU*C(t) exp ^F(t) + Ca (1 + TV(pi)) ^ rC,(r)C(r) dr^j Ufc - p 2 || L1 
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+tC(t) (l+C A j o TC^(T)C(T)n(T,t)dT^ 

X - f 2 )\\lj<x>([0,R]xM N *xKt) + \\^ 1Vx (f 1 ~ ^ 2 )llL^R w *)xL°°([0,.R]xK t ), 

completing the proof. □ 

We need below the following consequence of Kruzkov Theorem [8j Theorem 5] . 
Proposition 4.2 Let N x G N and T > 0. Consider the conservation law 



d t p + div^ f(t, x, p) = 
p(t,0) = p 



(4.4) 



with f G C°([0,T] x R N * x R;R N *); d p f , d p V x f and V%f continuous wherever defined;, 
dpi, diva; / G L°°([0, T] x M Nx x [-H,H]) for all H > 0. ylsstime taai p G (L^L^)^;!) 
is sttc/i £/ia£ = /or a.e. cc G R " \ B^N,c(0,d) for a given d > 0. Moreover, f(t,x,0) = 
/or a// i G [0,T] and x G R^. CaH p the Kruzkov solution to (fO|) and let K = 
su Pte[o,T] Then > f or allt G [°> T ]> = ® for a.e. x G R^ \-B RJV:c (0, d+Vt), 



where V — ||^p/|| Ij o ([o ) T]xR^x[-ii',X])' 

Above, / is assumed to satisfy the usual Kruzkov conditions, see [TTJ (HI)], or [5l [8] . The 
proof essentially relies on [8j Theorem 1]. 

Proof of Proposition Choose an x G R Nx \ B RNx (0, d + Vt). Let 5 > be such that 
fijjivz (x, <5) n B RNx (0, (2 + Vt) = 0, so that B RNx (x, 5 + Vt)f] B r n x (0, d) = 0. Applying [8, 
Theorem 1], with u = p and v = 0, we have that 



/ |p(t,x)|dx< / |p(x)|dx 

JB„ N „(x,5) JB mN „(x,S+Vt) 



'B rNx (x,S) JB RNx (x,S+Vt) 

hence pit) vanishes a.e. outside B r n x (0, d + Vt). □ 

Proof of Corollary \2.5l Fix any positive T G I. Let d be such that p vanishes outside 
B RNx (0,d) and call K = B RNx (0,d + VT). Let X 6 C~(M, [0, 1]) be such that x(x) = 1 for 
all x G /C. Define the convolution in the space variable /* = f * x Xi so that /* has compact 
support in x. Then, thanks also to the a priori bound (|2.9p . /* satisfies (f) on the interval 
[0, T]. Hence to the problem 

d t p + div x /* (t, x, p, p(t)) = 
p = tp(t,p, (Ap(i))(p)) 
p{0,x) = p(x) 
p(0) =p 

Theorem 12.21 can be applied, yielding the existence and uniqueness of a solution (p,p) in the 
sense of Definition 12.11 defined on all the interval [0,T]. Let now f(t,x,p) = f* (t, x, p,p(t)) . 
Then, p is a Kruzkov solution to ()4.4p and by Proposition 14.21 its support is contained in JC, 
for all t G [0,T]. Thereof ore, on the same time interval, by the definition of /*, (p,p) is the 
unique solution also to (jl.ip , always according to Definition 12.11 The rest of the proof easily 
follows. □ 
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